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ABSTRACT 



Context. PG 1535-1-547 is a bright Narrow Line Seyfert 1 galaxy, whose high energy emission shows strong variability both in shape 
■ and flux. On the basis of ROSAT observations, it is classified as "soft X-ray weak QSO", a class of objects whose X-ray-to-optical 

flux ratio is smaller than in typical QSOs. Their X-ray spectra are often characterized by highly ionized, complex absorbers and/or 
reflection from the inner accretion disk, and the relative importance of the two is currently debated. Whatever the correct interpretation 
may be, the presence of such features implies that we are looking at matter located in the innermost regions of these AGN. 
Aims. We want to clarify the nature of the X-ray emission of PG 1535+547, and constrain the physical properties of its innermost 



^ ' regions, where this emission originates. 

' Methods. We present new XMM-A'ewton observations of PG 1535+547 (90 ksec exposure time), from which we obtained two spectra 
separated by about one week, that we compare with a previous (about three years) XMM-Newton observation. 
Results. These observations support the complex and variable nature of the X-ray emission of PG 1535+547. The broad band observed 

^ ■ flux increases by a factor ~ 2.3 in three years, and then decreases by a factor ~ 1.3 in about one week. In the new EPIC spectra strong 

' absorption features at £ < 3 keV and a complex spectral shape in the iron line energy range are evident, coupled with a drop in 

I the emission at higher energies. We describe all the different states in a consistent way, assuming either a warm absorber plus a 

, relativistically blurred ionized reflection, or a two-phase warm absorber partially covering the source with the addition of a scattered 

■ component. 

• ' Conclusions. The observed variability can be ascribed mostly to warm absorbing gas in the innermost regions of PG 1535+547, that 

viS I appears to vary in its physical properties on timescales of both years and days. In the blurred reflection scenario all the analysed states 

^— ^ , require a high fraction of reflection from the disk, calling for some mechanisms able to increase the reflection component with respect 

25 ' to the intrinsic continuum. Finally, the strong variability observed in the X-ray band opposed to a more constant emission at optical 

' frequencies changes the value of the X-ray-to-optical spectral index, implying that PG 1535+547 can not actually be classified as a 

^ I soft X-ray weak AGN. 
■ i~H , 

^ Key words, galaxies: active - quasars: individual: PG 1535+547 - X-rays: galaxies 

1 . Introduction The inner region of AGN houses also "warm" (i.e., partially 

photoionized, T ~ 10^ K) gas that absorbs the nuclear emission 
Studying the X-ray emission from Active Galactic Nuclei in the X-ray band. This "warm absorber" (Halpern 1984) is of- 
(AGN) offers the potential of investigate the central engine pow- ten found to be complex and multi-phase, and outflowing with 
ering these sources. The fluorescent Fe emission line complex velocities in the range 0.001 - 0.01 c (e.g., Piconcefli et al. 2004; 
at 6.4 - 7 keV is an important direct probe of the dense mat- Ashton et al. 2004; Blustin et al. 2005; McKernan et al. 2007). 
ter around the nuclear region, from the molecular torus and the Recendy, more extreme outflow velocities (i.e., up to ~ 0.3 c) 
Broad Line Region (BLR) to the inner accretion disk, down to have been detected in the X-ray spectra of a few dozen of AGN 
a few gravitational radii (r^ = GMrh/c^)- In the latter case, in- (foj. ^ review, see Cappi 2006). The knowledge of the properties 
trinsically narrow emission lines emitted by the accretion disk of ^arm absorbers such as their location, covering fraction, mass 
are predicted to be broadened by Doppler and relativistic effects, content and outflow velocity can help us to understand not only 
Detailed studies of their energy profiles and flux variability could the accretion/ejection mode of AGN, but also their impact on 
give strong constraints on the geometry and kinematics of the ac- the suiTounding medium. Constraining the physical properties 
cretion disk in terms of inclination, inner and outer radii of the of ionized gas harbored in the inner regions of AGN and dis- 
emitting region and ionization state (e.g.. Miller et al. 2006; see entangle reflection from absorption contribution are challenges 
Reynolds & Nowak 2003 for a review). that can be unfolded only via high energy, temporaUy resolved 
specttal analyses. 

Send offprint requests to: L. Ballo, e-mail: Linked to both these topics, the origin of the so-called "soft 

lucia.balloOsciops.esa.int excess", an enhancement of flux seen below 2-3 keV with 
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respect to the extrapolation at low energies of the hard nuclear 
component, is still a debated point. While the constancy of its 
"temperature" over a wide range of luminosities and black hole 
masses seems in contrast with the idea of direct thermal emission 
from the accretion disk (see e.g. Piconcelli et al. 2005), two dif- 
ferent explanations have been proposed (both probably playing a 
fundamental role in a complex combination, see e.g. Chevallier 
et al. 2006): photoionized reflection from the inner accretion disk 
(Crummy et al. 2006), or absorption from a relativistically out- 
flowing, moderately ionized medium (Gierlinski & Done 2004, 
2006; Schurch & Done 2006; Sobolewska & Done 2007). 

Over the past years, a comparison between optical and X- 
ray properties of bright sources pointed out the existence of 
a population of AGN notably faint in soft X-ray relative to 
their optical fluxes, the so-called "soft X-ray weak QSO" (Laor 
et al. 1997). Analysing the X-ray emission of the Boroson & 
Green (1992) sample of nearby AGN, Brandt et al. (2000) found 
that approximately 10% of optically selected AGN are X-ray 
weak at soft energies, with an X-ray-to-optical spectral index 

ttox = log(f2kev/-F3oooA)/log(^2kev/v3oooA) < "2. The Strong 
correlation observed between qtox and absorption features in the 
UV band (mainly the C iv absorption equivalent width) supports 
the idea of X-ray absorption as primary cause of their soft X-ray 
weakness. This makes the X-ray weak AGN important target to 
study the ionization structure and the matter distribution of the 
warm absorbers, whose presence seems to be linked to absorp- 
tion in the UV band (Brinkmann et al. 2004; Piconcelli et al. 
2004). In this case, a sample composed by soft X-ray weak AGN 
would have a high incidence of X-ray obscured sources, playing 
an important role e.g. in the X-ray background synthesis models 
as well as in unified AGN schemes. 

On the other hand, Risaliti et al. (2003) from a "mini-survey" 
of 18 quasars observed with Chandra found a high incidence of 
intrinsicaUy under- luminous objects, suggesting that they rep- 
resent a population with a Spectral Energy Distribution (SED) 
difl'erent from that of standard blue quasars. So, a difl'erent in- 
trinsic X-ray emission mechanism cannot be completely ruled 
out. A strong variabiUty of the intrinsic continuum might be an- 
other reason for a X-ray weakness at the time of the observa- 
tion (see e.g. flie case of PG 0844+349; Brinkmann et al. 2003). 
Recently, Schartel et al. (2007) explained the weak X-ray emis- 
sion of PG 2112+059 with light bending (Miniutti & Fabian 
2004) instead of absorption, suggesting these sources as ideal 
target to observe the accretion disk near the central black hole. 

In this paper we present new XMM-Newton data of a soft X- 
ray weak QSO, PG 1535+547, comparing the results with pre- 
vious XMM-Newton observations. In our analysis, we assume a 
cosmology with Qm= 0.3, £1^= 0.7, and Hq=1Q km s~' Mpc"^ 

2. PG 1535+547 

PG 1535+547 (z = 0.038, De Vaucouleurs et al. 1991) is a 
radio-quiet (ratio of radio-to-optical flux density' R - 0.14; 
Kellermann et al. 1989) type I AGN comprised in the Palomar- 
Green Bright Quasar Survey (PG BQS; Schmidt & Green 1983, 
Mb = -20.8 for the cosmology adopted in the present paper). 
Because of the narrow Hy8 Hne (FWHM(Hy8) = 1480 km s ', 
Boroson & Green 1992; Veron-Cetty et al. 2001), the source is 
classified as Narrow Line Seyfert 1 (NLSl, which are histori- 
cally defined to have FWHM(I^) < 2000 km s"\ Osterbrock & 
Pogge 1987, and ratio of [Oiii]/!^ < 3, Shuder & Osterbrock 
1981). Its optical spectrum presents strong Feii emission lines 

' Measured at 6 cm and 4400 A, respectively. 



(de Veny & Lynds 1969; Phillips 1978; Smifli et al. 1997), as 
typically observed in NLSl galaxies (Boroson & Green 1992; 
Sulentic et al. 2000; Grape et al. 2004b). The optical/UV contin- 
uum appears substantially reddened (Smith et al. 1997) and po- 
larized: its optical polarization, P - 2.5%, is the highest among 
the PG sample (showing a {P) - 0.5%; Berriman et al. 1990). 
Strong intrinsic Civ /1 1549 absorption features, blueshifted by 
~ 2500 km s"' with respect to the redshift of the source 
(Sulentic et al. 2006) further classify the source as a "mini-Broad 
Absorption Line (BAL) QSO"^, and imply the presence of a 
moderate velocity outflow. Finally, no optical micro variability 
has been observed (Carini et al. 2007). 

The mass of the central black hole of PG 1535+547 has been 
derived by Vestergaard & Peterson (2006) and Zhang & Wang 
(2006) using two sUghtly difl'erent versions of the relation with 
the optical continuum luminosity and the FWHM(I^), and by 
Zhang & Wang (2006) by rescaling the stellar velocity disper- 
sion estimated from the [Om] line. The three evaluations give 
log(MBH/M0) = 7.19, 6.94 and 7.34, respectively. The slightly 
different Mbh derived from the FWHM(^S) and using the [O m] 
line width could suggests that the stellar velocity dispersion is 
larger than expected from the FWHM(H/S) - Mbh relation. This 
would imply for PG 1535+547 a location below the Mbh - o" 
relation found for local galaxies {e.g., Tremaine et al. 2002). 
A similar result was found for other NLSl galaxies accreting 
close to the Eddington hmit (Grape & Mathur 2004; Mathur & 
Grupe 2005), as in the case of PG 1535+547: the bolometric 
luminosity and the black hole mass estimated using the [O in] 
line from Zhang & Wang (2006) lead to an accretion rate rela- 
tive to the Eddington rate of logm = -0.32. It is worth noting 
that both NLSls and BAL QSOs are found to lie at the extreme 
of the anti-correlation between [O iii] and Fe ii emission (the so- 
called Eigenvector 1; Boroson & Green 1992) observed in low- 
redshift QSO, that is thought to be driven by a high Eddington 
ratio (Boroson 2002). 

In the X-ray band, ROSAT PSPC observations performed on 
1993 detected the source with a 0.1 - 0.4 keV count rate^ of 
(0.761 ± 0.195) X 10"^ cts s"'. The source was not detected in 
the hard (0.5-2.0 keV) band, nor it is detected in the ROSAT AW 
Sky Survey catalog^. Gallagher et al. (2001) estimated an upper 
limit to the observed count rate in the 0.5 - 2 keV energy range^ 
of 3.2 X 10"^ cts s"'. The non-detection above 0.5 keV reminds 
anoflier NLSl, the X-ray transient WPVS 007. In their analy- 
sis of ROSAT data, Grupe et al. (1995) suggested obscuration as 
the possible cause of the turn-off observed in the 0.5 - 2.4 keV 
energy range. On the basis of the ROSAT observation, Brandt 
et al. (2000) classified PG 1535+547 as soft X-ray weak QSO. 
In 1999, six years after the ROSAT observations, ASCA found 
that the source emission has varied in a significant way: assum- 
ing the best fitting ASCA spectral model, Gallagher et al. (2001) 
predicted a PSPC count rate of 4:9^^^ x 10"^ cts s"', i.e. a fac- 
tor ~ 1 .5 higher than the upper limit estimated for the ROSAT 
observation. Moreover, the authors found evidence of X-ray in- 
trinsic absorption (A^h ~ lO^'^ cm"^) with only partial covering 



^ Defined as QSO showing in their optical/UV spectrum absorption 
lines with intermediate widths between the "classical" BAL QSO (ab- 
sorption line width > 2000 km s"') and the Narrow Absorption Line 
QSO (absorption line width < 500 km s '); see e.g. Weymann et al. 
(1991); Narayanan et al. (2004). 

' ftp://ftp.xray.mpe.mpg.de/rosat/catalogues/lrxp/ 
http://www.xray.mpe.mpg.de/cgi-bin/rosat/src-browser 

' The authors reported the detection of 5 photons above 0.5 keV in a 
1.'2 radius source cell, with 2 expected from the background. 
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Table 1 EPIC XMM-Newton data observation log. 



Obs.ID 


Instrument 


Filter 


Date 


Start 
End 


Exp. 
[ksec] 


Net Exp. Time 
[ksec] 


Net Count Rates" 
[10"' counts sec"'] 


S/N" 


Name used 
in the text 




MOSl 


thin 


2006 Jan 24-25 


22:16:06 
05:26:53 


25.7 


20.64 


9.5 ±0.2 
0.5 ±0.1 


43.94 
10.20 




0300310501 


M0S2 


thick 


2006 Jan 24-25 


22:16:06 
05:26:58 


25.7 


21.60 


8.6 ± 0.2 
0.5 ±0.1 


42.68 
9.87 






pn 


thin 


2006 Jan 24-25 


22:38:48 
05:27:13 


24.0 


15.61 


34.1 ±0.5 
2.7 ±0.1 


72.01 
20.04 


(A) 




MOSl 


thin 


2006 Jan 22-23 


22:43:11 
06:53:58 


29.3 


18.31 


9.9 ± 0.2 
0.5 ±0.1 


42.11 
9.88 




0300310401 


M0S2 


thick 


2006 Jan 22-23 


22:43:11 
06:54:03 


29.3 


18.42 


8.4 ± 0.2 
0.5 ±0.1 


39.04 
3.91 






pn 


thin 


2006 Jan 22-23 


23:05:53 
06:54:18 


27.6 


13.77 


30.4 ± 0.5 
2.7 ±0.1 


63.82 
5.60 






MOSl 


thin 


2006 Jan 16-17 


22:36:41 
07:17:28 


31.1 


8.38 


14.1 ±0.4 
0.4 ±0.1 


30.15 
5.56 




0300310301 


M0S2 


thick 


2006 Jan 16-17 


22:36:41 
07:17:33 


31.1 


12.93 


12.3 ± 0.3 
0.6 ±0.1 


39.54 
8.84 


(B) 




pn 


thin 


2006 Jan 16-17 


22:59:23 
07:17:48 


29.4 


5.80 


46.7 ± 0.9 
3.2 ± 0.2 


51.79 
13.48 






MOSl 


thin 


2002 Nov 03 


01:06:17 
09:20:24 


29.5 


24.53 


4.0 ±0.1 
0.3 ±0.1 


30.71 
7.98 




0150610301 


M0S2 


thin 


2002 Nov 03 


01:06:17 
09:20:24 


29.5 


24.53 


4.1 ±0.1 
0.3 ±0.1 


30.76 
7.87 


Nov 2002 




pn 


thin 


2002 Nov 03 


01:28:48 


27.8 


20.66 


13.6 ±0.3 


49.58 










09:20:44 






1.8 ±0.1 


18.16 





" In the energy range 0.4 - 10 keV (first line) and 6.5 - 10 keV (second line). 



of the power law continuum. The limited quality of the ROSAT 
data prevented the authors from ascribing the observed variation 
to an increase of the continuum flux or a decrease of the absorp- 
tion. As noted by Gallo (2006), NLSl galaxies with complex 
high-energy spectra (such as PG 1535-1-547) appear to be X-ray 
weak AGN. 

The XMM-A'ewton observation performed in 
November 2002 found PG 1535-H547 in a higher state than 
previous observations at energies E > 2 keV (Schartel et al. 
2005). Moreover, the EPIC spectra clearly require the presence 
of complex absorbing structure. A good description of the data 
is obtained only assuming a combination of an ionized absorber 
and a neutral absorber with a large covering factor. Finally, the 
data suggest the presence of a broad emission feature, whose 
profile is well described by a relativistically broadened Fe Ka 
emitted in a disk accreting around a black hole. The absence 
of a features like this in the ASCA data allowed the authors to 
suggest a possible variability in the emission Une. A similar 
behaviour (i.e., complex X-ray emission and strong variability) 
has been reported for other NLSls: for example, Grupe et al. 
(2007) found Mrk 335 in an extremely low X-ray state, whose 
interpretation includes partial-covering absorption or X-ray 
reflection from the disk. A comparison of the emission observed 
in 2007 with all previous observations shows a decrease in flux 
by a factor of more than 30, that the authors ascribe mainly to a 
change of the intrinsic absorption. 

3. Observation and data reduction 

We observed PG 1535-1-547 (RA=15:36:38.3, Dec= -1-54:33:33) 
with XMM-Newton in January 2006 in three pointings, the first 
separated by about one week from the others, for a total of about 
90 ksec (Obs. ID 0300310301, 0300310401 and 0300310501). 
The observations were performed with the European Photon 
Imaging Camera (EPIC), the Optical Monitor (OM) and the 
Reflection Grating Spectrometer (RGS; den Herder et al. 2001); 
the source is not detected with the last instrument. The three 



EPIC cameras (pn, MOSl, and MOS2; Striider et al. 2001; 
Turner et al. 2001) were operating in full frame mode, MOSl 
and pn with the thin filter and MOS2 with the thick filter apphed. 
The observation details are reported in Table 1. OM (Mason 
et al. 2001) performed 6 observations for each pointing, three 
wifli the UVM2 (Ae/f = 2310 A) and three wifli the UVW2 
(Aeff = 2120 A) filters in the optical light path, of 1 .9 ksec each; 
all OM exposures were performed in the "Science User Defined" 
image mode (with a windows size of 5' x5')- The XMM-Newton 
data have been processed using the Science Analysis Software 
(SAS version 7.0) with the calibration from February 2007; the 
tasks epproc and emproc were run to produce caUbrated and con- 
catenated event lists for the EPIC cameras. 

EPIC event files have been filtered for high-background time 
intervals, following the standard method consisting in rejecting 
periods of high count rate at energies > 10 keV. Events cor- 
responding to patterns 0-12 (MOS1&2) and 0-4 (pn) have 
been used (see the XMM-Newton Users' Handbook; Ehle et al. 
2001). We have also generated the spectral response matrices at 
the source position using the SAS tasks arfgen and rmfgen. We 
can exclude event pile-up in the MOS and pn data. Source counts 
were extracted from a circular region of radius 18" and 21" 
for the MOS and pn. Background counts were extracted from 
a nearby source-free circular region of 36" and 42" radii, re- 
spectively. The net exposure times at the source position, the net 
count rates and the S /N ratios in the energy ranges 0.4-10 keV 
and 6.5-10 keV are reported in Table 1 . 

We tested short time variability within each observation, 
generating source light curves in diff'erent energy ranges with a 
binning time of 5 ksec. No bin shows significant deviation from 
the mean value. 

We processed OM data with the omichain routine of SAS; 
for each OM filter and for each pointing three images were ob- 
tained. The count rate informations were extracted from the three 
omichain output and averaged, assuming the standard deviation 
as an estimate of their error. The count rates were converted to 
magnitudes (in the AB system) and fluxes using the conversion 
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Table 2 OM integrated magnitudes and fluxes (January 2006 observa- 
tions). 



Obs. ID 



AB magnitudes 
[mag] 



UVM2 
[2310 A] 



UVW2 
[2120 A] 



Flux 

[10"'^ ergs cm~^ s~ 



A '] 



UVM2 
[2310 A] 



UVW2 
[2120 A] 



0300310501 16.567 ±0.016 16.719 ±0.027 4.811 : 
0300310401 16.510 ±0.015 16.648 ±0.026 5.072: 
0300310301 16.496 ± 0.015 16.625 ± 0.025 5.139: 



0.069 4.977 ±0.123 
0.071 5.312 ±0.127 
0.071 5.427 ±0.127 



factors reported by the SAS homepage*'; the errors have been 
calculated propagating the uncertainties on the count rates. The 
obtained values are provided in Table 2; in all the observations, 
the ratio of the maximum to the minimum count rates is consis- 
tent with unity within the errors. 

In order to better compare the old and new XMM-Newton 
observations, we have re-analysed the EPIC XMM-Newton data 
taken in November 2002 (archived under the observation identi- 
fier 0150610301). At this occasion, PG 1535H-547 has not been 
observed with OM. Following Schartel et al. (2005), for the two 
MOS exposures the background was determined from an annu- 
lus around the extraction area of the source counts (inner and 
outer radii: 48" and 84", respectively). Details of this observa- 
tion, as well as exposure times, count rates and S /N ratios in the 
energy ranges 0.4 - 10 keV and 6.5 - 10 keV, are reported in 
Table 1. 



4. The X-ray emission: a first comparison between 
the datasets 

We performed the spectral analysis of the MOS and pn spec- 
tra over the energy range from 0.4 keV to 10.0 keV. The EPIC 
spectra have been analysed using standard software packages 
(FTOOLS version 6.2, XSPEC version 12.3. Iq, XSPEC ver- 
sion 11.3.2ag; Arnaud 1996). All the models discussed in the 
following assume Galactic absorption with a column density of 
A^H.Gd = 1.35 X lO^" cm-2 (Murphy et al. 1996). Unless other- 
wise stated, the figures are in the rest-frame of PG 1535-1-547, 
and fit parameters are quoted in the same frame. 

Comparing the three observations of January 2006 it is ev- 
ident that the flux decreases considerably passing from the ob- 
servation 0300310301 to the other two, that present the same 
spectral shape and intensity. This behaviour is particularly en- 
hanced at energies lower than 0.7 keV. Therefore we combined 
the pn and MOS data from 0300310401 and 0300310501 (the 
label A identifies the combined spectra), keeping separate the 
0300310301 data, that are labeled as B. In order to apply the;^'^ 
statistics, source counts for the (A) and (B) spectra were binned 
to have at least 40 and 30 counts in each energy bin, respectively. 

From a comparison between the new and the old datasets, 
a strong change of both flux and spectral shape is evident (see 
Fig. 1, left panel). 

In Fig. 1, right panel, we have plotted the ratio of the 
January 2006 XMM-Newton spectra of PG 1535+547 to the 
November 2002 data (considering that the MOS spectra agree, 
within the errors, with the EPIC-pn spectra, to avoid clutter only 
the last are shown). The best fit found by Schartel et al. (2005) 
for the November 2002 data (an intrinsic power law with super- 
imposed a relativistically broadened line emitted by an accretion 
disk around a Kerr black hole, covered by a combination of an 




^ http://xmm.vilspa.esa.es/sas/7.0.0/watchout/ 
Evergreen_tips_and_tricks/uvflux.shtml 



Fig. 1 Left panel: EPIC-pn spectra of the three observations: red filled 
circles, spectrum (A); grey open triangles, spectrum (B); sky-blue stars, 
November 2002 data. Right panel: Ratios between the new EPIC-pn 
spectra and the best fit model for the November 2002 dataset (Schartel 
et al. 2005). The black dotted line corresponds to the position of the 
6.4 keV fine. 



ionized absorber and a neutral absorber with a large covering 
factor) cannot describe neither the (A) data or the (B) spectra. 
This figure may offer a first indication about the origin of the 
observed emission. The strongest variation is observed in the 
soft energy band, with the emission between 0.7 and 1.5 keV 
being less affected by changes. This ratio would imply that the 
observed variability is associated mainly with absorbing matter, 
whose physical properties are different among the observations. 



5. Spectral analysis of the January 2006 (A) data 

We started with the (A) data, due to the higher S IN of these 
spectra. The analysis of the November 2002 and the (B) datasets 
was performed taking into account the results of this first part of 
the work. 



5.1. The emission at energies > 3.5 keV 

First, we fitted the (A) data with a simple power law. Due to the 
complex spectral shape both at low and high energies, we con- 
sidered only the 3.5 -5 keV energy range, finding a photon index 
r = 1.40 + 0.19. This spectral index appears flatter than the mean 
value found at energies > 2 keV for a sample of AGN selected 
in the PG BQS ((F) = 1.87 + 0.10; PiconcelH et al. 2005). This 
discrepancy is even stronger considering only NLSl galaxies, 
whose photon indices appear to be on average steeper than those 
of Seyfert galaxies with broad optical lines (e.g., studying ASCA 
spectra of 23 NLSls, Leighly 1999 found a mean value in the 
2-10 keV energy range (F) = 2.19 + 0.10). When the model is 
extrapolated to the whole 0.4 - 10 keV band, the most evident 
feature is a strong decrease in the flux below 3 keV, particularly 
enhanced between 0.7 and 2 keV (Fig. 2, upper panel). 

Focusing at £■ > 3.5 keV, this simple model shows residual 
structures in the 5.5 - 6.5 keV rest frame region, and is also 
unable to describe the high energy end of the spectra (Fig. 2, 
lower panel). 

To model the observed residuals, as a first step we added one 
or more narrow or broad line(s) to the powerl-law continuum 
(ZGAUSS component in XSPEC, with cr = 1 eV or free to vary), 
and/or an absorption edge at ~ 7 keV (ZEDGE component in 
XSPEC). The addition of the absorption edge can't reproduce 
the data between 6 and 7 keV, and still required an unusually 
low photon index, F ~ 1.27. The observed spectra can be well 
described assuming a model with two Gaussian lines, one nar- 
row and one broad, whose energies were consistent with being 
produced by neutral iron, or a power law model with superim- 
posed a disk line in a KeiT configuration (LAOR component in 



L. Ballo et al.: X-ray variability in PG 1535+547 



5 



Table 3 Analysis of the high energy (3.5 - 10 keV) part of the January 2006 (A) pn and MOS spectra. 
Results for a reflection-based model. 

Power law+Gaussian^ Kdblur[Reflion] ^ 

Model r Norn? E^f [keV] ^ R~jGMjF] Tf] ^Vd.o.f. 

Fe/solar Norm^ 

disk+line 2.32:;;-^° 6.401; »^ 6.46 ± 0.04 3.02+" ^^ 3.85!^-!,° 0.02™{;] 109.2/119 

^•^-'-0.22 ^-^'^•-^'^-127.41 ^-"^"-0.17 

Results for an absorption-based model. 

Power law+Gaussian^ PC absorber Absorber 

Model r Norm" % [keV] ' A'h" Cov. fraction"" A^h" ;r'/d.o.f. 

neutral pc+line 2.55!;;^= 40.871j^;»^ 6.411^^^ 10.79!^;^^ OM^g 3.241!;\'» 121.3/123 

EiTors are quoted at the 90% confidence level for I parameter of interest (Ax^ = 2.71). ^ Narrow line width fixed to e V. Photon index coupled with the primary power law photon 
index; outer radius fixed to the maximum value, /?out ^ ^S' " units of 10""^ photons keV"^ cm"- s"' @ 1 keV. Ionization parameter ^ = 4;rFjjj /« (erg cm s"' ), where Fjjj is 
the ionizing flux and n is the hydrogen nucleus density (part cm""*) of the illuminated slab. ' In units of 10"^ photons keV"' cm"- s"' @ 1 keV of the reflected spectrum. ^ In units of 
10^^ cm"^. ^ Complement to one of the covered-to-uncovered flux ratio. 




Energy rest frame (KeV) 



Fig. 2 Data/model ratio for the (A) spectra (green open squares, MOS; 
red filled circles, pn). The adopted model is a simple power law with 
Galactic absorption fitted between 3.5 and 5 keV, and extrapolated up 
to 10 keV (lower panel), and to the whole studied energy range (0.4 - 
10 keV, upper panel; for demonstration purposes, here the data have 
been binned to have a significance of 8cr up to 30 bins). In both panels, 
the black dotted line marks the position of the 6.4 keV line. 

XSPEC; Laor 1991; Fabian et al. 2002). The addition of a nar- 
row Gaussian line with energy consistent with being produced 
by matter in low ionization stages [i.e., Fe i-Fe xvii) improves the 
fit significantly (F-test probability' 97.3%). Being the line unre- 
solved by EPIC-pn, in the following its width has been fixed to 
the instrumental spectral resolution. 

In the last years, two different scenarios have been proposed 
to account for high-energy spectral shapes as observed in the (A) 
spectra (Fig. 2, lower panel): a "reflection-based" (e.g., Mrk 335, 
Longinotti et al. 2007; O'Neifl et al. 2007; HE 0450-2958, Zhou 
et al. 2007) or an "absorption-based" (e.g., IH 0707-495, Boiler 
et al. 2002; Gallo et al. 2004; Fabian et al. 2004; Tanaka et al. 
2004; IH 0419-577, Pounds et al. 2004; Fabian et al. 2005; 
IRAS 13197-1627, Dadina & Cappi 2004; Miniutti et al. 2007) 
one. Despite their different shape in the Fe line energy region, 
the superposition of different layers of partially ionised absorb- 
ing material could mimic the presence of a broad line component 

' But see the caveats in using the F-test to measure the significance 
of narrow lines described in Protassov et al. (2002). 



superimposed to the continuum, a feature naturally present in 
relativistically blurred ionized reflection model (see e.g. Reeves 
et al. 2004; Grupe et al. 2007; Petrucci et al. 2007). 

A reflection component can reproduce in a good way 
(X^/d.o.f. - 109.2/119) the high-energy curvature only if orig- 
inates in the accretion disk**, leading to a rather steep intrinsic 
continuum, F ~ 2.3. To model the ionized reflection in XSPEC 
we used the REFLION"^ external table, which incorporates both 
line emission with Compton broadening and reflection contin- 
uum (Ross & Fabian 2005; Ross et al. 1999). Then, the reflec- 
tion from the disk has been constructed convolving this compo- 
nent with the relativistic blurring model KDBLUR, produced by 
A. Fabian and R. Johnstone using the same kernel of the LAOR 
disk line model'". 

The 3.5 - 10 keV data can be equally well reproduced 
(X^/d.o.f. - 121.3/123) in an absorption scenario, assuming a 
high column density (A^h ~ 10^^ cm"^) cold absorber covering 
~ 80% of the primary continuum. This latter component is rather 
steep (F ~ 2.5), to account for the observed spectral drop at the 
highest energies (E > 1 keV). Allowing the ionization state of 
the absorbers or the iron abundance to be free parameters, the fit 
does not improve further. 

The high-energy curvature observed in the (A) spectra of 
PG 1535-1-547 can thus be well accounted for in both reflection- 
based and absorption-based scenarios (see Table 3 and Fig. 3). A 
gaussian emission line is detected with an F-test probability of 
98.6% in both cases; its energy and EW (~ 50 - 80 eV) are com- 
patible with Ka fluorescence emission of almost neutral iron, 
produced by transmission through cold gas of column density 
< 10^^ cm"^ totally covering the source (Makishima 1986). 

5.2. The whole energy range 

We then extended the analysis at all energies (0.4 - 10 keV). The 
(A) spectra between ~ 0.7 and 2 keV appear strongly absorbed, 
with an important component rising in the softer band. Evident 
absorption features are present at £ ~ 2 keV as well as at ener- 



Reflection from neutral material far from the nucleus (PEXRAV in 
XSPEC; Magdziarz & Zdziarski 1995) does not provide a good descrip- 
tion of the data CrVd.o.f. = 132.2/121). 
' http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/models/reflion.html 

Unless otherwise stated, the outer radius was fixed to its maximum 
value, i^out = 400 rg, while the inner radius and the emissivity power 
law index were free to vary. 
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Table 4 Analysis of the broad band January 2006 (A) pn and MOS spectra (0.4 - 10 keV). 

Result for a warm absorber* [pailial covering warm absorber*(intrinsic power law+disk reflection)+narrow line] model. 

Power law PC warm absorber Warm absorber Kdblur[Reflion] T 

Obs. r Norm" f Cov. fraction" N^' f ^8 Rj^ [fg] p ;f^/d.o.f 

Fe/solar Refl. fraction*" 

(A) 2.06±0.08 OF™ 6.11;:^f'3'^ 68.14+^^77 OW^^ 2.13±0.08 3.67±0.08 434^ 1.31+|!,-j;:^ 121.59+*'22 290.1/281 

0.003 ± 0.01 4.54+'j'-35 65.2% ' 

Result for a neutral absorber*{warm absorber* [partial covering warm absorber*(intrinsic power law)]+scattered component+narrow lines} model. 

Power law'f PC warm abs. Warm abs. Gussian Neutral abs. 

Obs. T Norm" f Cov. fraction" N^' f N^' ^^/d.o.f 

Tcran Scatt. fraction'' 



scatt 

(A) 2A3t''o'oi 2l74^H? 8.43+|5^ 19.96+^J^ 0.43™y-; 4.49™;^ 94.04+y-^» 039^^ <3() 03F™ 297.5/279 

T 00+039 n nnf,+8:002 



EiTors are quoted at the 90% confidence level for 1 parameter of interest (Ax^ - 2.71); narrow line energy fixed to E^^ - 6.4 keV, narrow line width fixed to eV. ' Photon index 
coupled with the primary power law photon index; outer radius fixed to the maximum value. i^Qjjf = 400 rg. First line, intrinsic power law; second line, scattered power law. 
" In units of 10"'* photons keV"^ cm"^ s"' @ 1 keV. ^ In units of 10"^ cm"". Ionization parameter f = 4;rF-jj/» (erg cm s"' ), where F^jj is the ionizing flux and n is the hydrogen 
nucleus density (pait cm"-^) of the illuminated slab. ^ Complement to one of the covered-to-uncovered flux ratio. ^ Ratio between the ionized reflection flux and the total flux in the 
(0. 1 — 1000) keV energy range. ^ Line centroid energy, in keV. ^ Line width, in eV. ^ Ratio between the scattered flux and the primary, unabsorbed flux. 
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Fig. 3 EPIC-pn January 2006 (A) spectra of PG 1535+547 in the 
3.5-10 keV energy range, plotted in comparison to the reflection-based 
(left panel) and the absorption-based (right panel) models (see §5. 1 and 
Table 3). 



gies lower than about 0.7 keV. Simple models, composed by a 
single power law continuum covered by a combination of neu- 
tral or ionized absorbers (modeled with ABSORI component in 
XSPEC; Zdziarski et al. 1995), were not able to reproduce the 
observed shape. As obtained considering only the high-energy 
part of the spectra, a model based on reflection from neutral ma- 
terial distant from the black hole alone was not able to explain 
the observed emission, too. More complex scenarios have to be 
investigated. 

The basic models adopted to test the reflection-based and the 
absorption-based scenarios are the ones described in §5.1: 

- a power law plus reflection from the disk (the latter con- 
structed in XSPEC using the REFLION external table con- 
volved with the KDBLUR model); 

- a simple power law and a complex partial covering material. 

The overall broad-band shape was roughly reproduced by the 
two models, although in the latter case the strong soft excess 
observed can be partially accounted for only allowing the ab- 
sorbing column densities to be ionized. Under this assumption, 
the absorbers are able to let the nuclear flux to be transmitted at 
low energies. However, these simple models were inadequate to 
model the absorption features, resulting in statistically unaccept- 
able fits. 



For the reflection-based model, following the best fit ob- 
tained in §5.1 we adopted a continuum model composed by an 
intrisic power law and a relativistically blurred ionized reflec- 
tion. To model the absorption features observed in the spectra, 
we assumed the presence of a warm absorber partially cover- 
ing the central region; a second distribution of matter with lower 
ionization covering the whole emission is also required. Finally, 
a transmitted narrow line at 6.4 keV (rest frame) completes the 
model. This complex model resulted to be a good description of 
the (A) data Cr-/d.o.f. - 290.1/281); the best fit parameters are 
shown in the first part of Table 4. 

As previously noted, the presence of ionized absorption can 
affect the high-energy spectrum and the iron line profile. We 
then checked our results modeling the absorber in a more so- 
phisticated way, adopting a more refined photoionization model, 
namely XSTAR' ' (Kallman & Bautista 2001). XSTAR is a pub- 
lic code that consistently (i.e., by solving temperature and ion- 
ization balance equations) calculates the opacities of a spheri- 
cal shell of Compton-thin gas photoionized by a central source 
of radiation. It offers the possibility to construct tables varying 
the photoionizing continuum spectral shape and strength, the gas 
physical conditions (density, pressure, thickness) and the abun- 
dances of cosmic elements. Its atomic database includes all the 
most relevant transitions. To perform the fit using XSPEC ver- 
sion 12'^, we considered A^h to be a free parameter but keeping 
its value fixed during the minimization. The best fitting param- 
eters are in good agreement with the values reported in the first 
part of Table 4, while the fit obtained with this change in the 
absorption model is statistically worse. Above all, a disk com- 
ponent is still statistically required. 

Regarding the two-phase warm absorber-based scenario, to 
improve the fit we added to the power law (covered by an ion- 
ized, partial covering absorber and an ionized, fully covering ab- 

' ' http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html 

The model assumed to describe the reflection from the disk taking 
into account the relativistic blurring can be fitted only inside the new 
version of XSPEC. On the other hand, using the XSTAR tables in the 
last version of XSPEC, the minimization algorithm suffers problems 
in interpolating the column density values between grid points. Their 
origin is currently under investigation (Craig Gordon, priv. comm.). 



L. Ballo et al.: X-ray variability in PG 1535+547 



7 



Table 5 Simultaneous analysis of the three available XMM-Newton datasets for PG 1535+547 in the 0.4 - 10 keV energy range: the adopted 
model was the accretion disk reflection model (a continuum composed by an intrinsic power law and a reflection from the disk partiafly covered 
by an ionized absorber; the whole continuum was seen tlirough a warm absorber; ;(;^/d.o.f. = 379.6/364). Luminosity are absorption-corrected, 
while fluxes are corrected only for the Galactic absorption. 



Power law PC warm absorber Warm absorber Kdblur[Reflion] ^ 



Obs. 


r 


Norm" 






Cov. fraction" 






in 


Fe/solar 


Refl. fraction" 


(A) 


2 05+""' 

'^■"•'-0.12 


2 76+" '*'* 


5 77+'-'*'' 

J.I '_i9Q 


5C 1Q+25.S6 


0.72 ± 0.08 


^■"'-0.28 


n 9O+0.70 

3 -^S- 1.46 


4 26+"-" 

^■^°-0.27 

(2.7 ± 0.6) X 10-3 


1.31+1-^" 

4 7i+8:§? 

^•'^-1,83 


130 71+'"" 
^■'"■'^-19.31 

48.8% 


(B) 


2.05* 


T C7+0.55 
^••"-1.01 


a 44+2.13 


66.62+g:BB 


0.85 ± 0.06 


2.07* 


3.28* 


4.26* 
(2.7 X 10-3)* 


1.31* 
4.71* 


125 97+12.14 
60.1% 


Nov. 2002 


1 7Q+0.10 
■' -0.09 


1 55+0-35 

-0.96 


20.94:^;« 


98.87:30.08 


-0.01 


"■■'^-0.36 


3.28* 


4.26* 
(2.7 X 10-3)* 


1.31* 
4.71* 


65 43+21.50 

"•'■^■'-17.23 

45.0% 



Observed flux [10 '3 ergs cm ^ s '] Unabsorbed luminosity [10^3 ergs s '] 

Obs. 0.4 - 10 keV 0.4 - 1 keV 1 - 3 keV 3 - 10 keV 0.4 - 10 keV 0.4 - 1 keV 1-3 keV 3-10 keV 

(A) 25.54^;"" 1.52++^+ 3.11+;-^* 20.90!^'^^^ lOT 088 039 080 

(B) 33.53!^|| 2.74:f30 4.16:!:^-i| 26.63:^21.83 2.6I 1.18 0.44 1.00 
Nov. 2002 15.U+\°^* 0-37;^;f^ ^-^^In^l 13.84;::^;°^ 0.96 014 019 0.63 



Errors are quoted at the 90% confidence level for 1 parameter of interest (Ax = 2.71); narrow line energy tixed to E^^ = 6.4 keV. For the accretion disk reflection component, 
the outer radius was fixed to the maximum value, /?out ~ 400 rg. ^ For the single state, the power law and the accretion disk reflection components have the same photon index. 
^ Parameters for the states tied in the fit. " In units of 10"" photons keV"' cm-^ s-' @1 keV. In units of 10^^ cm-^. ' Ionization parameter f = 4;rFjjj/n (erg cm s-'), where Fyj 
is the ionizing flux and n is the hydrogen nucleus density (part cm'') of the illuminated slab. ^ In units of rg. Complement to one of the covered-to-uncovered flux ratio. ' Ratio 
between the ionized reflection flux and the total flux in the (0.1 - 1000) keV energy range. 



sorber) a soft power law component absorbed only by neutral 
gas (A;^'^/Ad.o.f. - 14.9/1). If we let the soft power law photon 
index to be a free parameter, the fit improves significantly but 
still remains statistically poor (x^/d.o.f. = 350.0/279). Looking 
at the residuals, it seems that the ABSORI model is inadequate 
to describe the complex shape of the absorbing structure seen at 
0.7 - 1 keV. Moved by these considerations, we try to reproduce 
the soft complexities with the cited XSTAR tables. 

In order to apply these tables, for this part of the analysis we 
passed to XSPEC version 1 1 . We fitted to the data our model 
(a double ionized absorber, one partially covering and the other 
totally covering the source, plus a soft scattered component) us- 
ing a pre-calculated XSTAR table. The publicly available table 
"grid 18b" assumes a 10'*'* ergs s"*, F = 2 power law ionizing 
a constant density shell of gas whose abundances can be varied. 
The turbulent velocity of gas is assumed to be 100 km s"*. We 
obtain a good fit to the data with an oxygen overabundance of 
~ 4xsolar, requested to reproduce the deep absorption structure 
observed between 0.7 and 1 keV. We note here that an ionizing 
continuum steeper than F = 2 results in deeper resonant absorp- 
tion in this energy range due to the higher number of soft ion- 
izing photons (see e.g. Nicastro et al. 1999). If a grid generated 
assuming a continuum with F - 2.5 is applied to the model, no 
oxygen overabundance is needed to reproduce the depth of the 
absorption structure; however, in this case the enhanced opac- 
ities at low energies {E < 0.7 keV) caimot reproduce the flat 
spectra observed in this energy range. We decided to proceed 
with the analysis adopting the "grid 18b" table, with these con- 
siderations in mind. 

A visual inspection of data/model ratio shows an excess of 
emission in the soft band. The addition of a gaussian emis- 
sion line improves the fit. The line energy, E ~ 0.59 keV, is 
fully compatible with the photoionized O vii triplet energy. Soft 
X-ray emission lines, mainly Ovii or Neix, have been found 
in the spectrum of other Seyfert galaxies (e.g., Mrk 6, Immler 
et al. 2003; Mrk 1239, Grupe et al. 2004a; Mrk 304, Piconcefli 
et al. 2004; Guainazzi & Bianchi 2007 and references therein). 



Broad components of the O vii triplets associated with the BLR 
have been observed in the X-ray spectra of the bright NLSl 
Mrk 1 10 (Boiler et al. 2007) and of the Seyfert galaxy Mrk 279 
(Costantini et al. 2007). Unfortunately, the flux of PG 1535+547 
is too low and prevents us from drawing information from the 
RGS spectra. 

The best fit model results in a moderately ionized gas (A^h ~ 
8.4 X 10^2 cm ^ ~ 19.9 erg cm s"*) covering ~ 40% of 
the primary emission, a second absorber with lower column 
density and higher ionization state (A^h ~ 4.5 x 10^^ cm"^, 
^ ~ 94 erg cm s ') fully covering the source, a scattered com- 
ponent and two gaussian emission lines (Fe and O vii). The ab- 
sorbed power law is rather steep (F - 2.5), and the scattered one 
is even steeper, as expected if scattering occurs in a photoionized 
gas. Finally, the whole emission is covered by a neutral absorber 
with a rather low column density (A'h ~ 2 x 10^* cm"^). This 
complex model results in an equally good description of the (A) 
spectra (x^/dof- 297.5/279, see second part of Table 4). For 
comparison, in Fig. 4 we have plotted the unfolded best fits to 
the pn data obtained with a partially absorbed accretion disk re- 
flection (left panel), or a combination of ionized absorbing layers 
and the addition of scattered emission (right panel). 

A partially ionized absorbing material with large velocity 
shear (the SWIND'^ model in XSPEC, originally proposed by 
Gierlihski & Done 2004 and updated by Gierlihski & Done 
2006) as origin of the observed soft excess cannot be tested here, 
being the maximum column density assumed in the analytical 
model (A^H ~ 8 x 10^^ cm"^) lower than the value we need to 
reproduce the high energy spectral shape. 

6. The origin of X-ray variability: comparison 
among the different XMM-Neivton datasets 

The physical scenarios explored before to describe the (A) data 
yielded two complex models, composed by several continuum 

'3 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/models/swindl.html 
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Table 6 Simultaneous analysis of the three available XMM-Newton datasets for PG 1535+547 in the 0.4 - 10 keV energy range: the adopted 
model was the complex absorption scenario: an ionized, partially covering absorber plus a scattered component emerging from an ionized 
layer that fully covers the source; ;if^/d.o.f. = 417.5/376). Luminosity are absorption-corrected, while fluxes are corrected only for the Galactic 
absorption. 



Power law^^ PC warm abs. Warm abs" Gaussian Neutral abs. 



Obs. 


r 


Norm" 






Cov. fraction" 














'"scatt 


Scatt. fraction* 
















(A) 
(B) 


^■■'^-0.49 
2 90+0.65 

2.54* 
2.90* 


32.56+^^,f 
005+"-'""' 

"■""•'-0.002 

37.33:!:5;^5 

0.004 ± o!ooi 


8.24^^:" 
8.24* 


18.16!^:^^ 

18.16* 


0.62 ± 0.03 
51+°-'^ 

-0.03 


a QC+0.42 
•^■'•^-0.23 

3.95* 


91 32+'^-'*' 
91.32* 


0.56 ± 0.02 40 ± 30 
0.56* 40* 


"■^"^-0.01 

0.22* 


Nov. 2002 


2.54* 
2.90* 


30.60!f3« 

"■""'-0.003 


8.24* 


18.16* 


91+0-03 

"•'^-0.04 


3.95* 


60 69+i''-io 


0.56* 


40* 


0.22* 



Observed flux [10 " ergs cm " s '] Unabsorbed luminosity [10^^ ergs s '] 

Obs. 0.4 - 10 keV 0.4 - 1 keV 1 - 3 keV 3 - 10 keV 0.4 - 10 keV 0.4 - 1 keV 1 - 3 keV 3 - 10 keV" 

(A) 25.30"^" 1.53 ±0.01 3.18;"-'» 20.59"-^^ 5?70 172 TST Ul 

(B) 33.15*™ 2.60 + 0.01 4.48+|* 26.06+f| 7.24 3.55 2.31 1.38 
Nov. 2002 14.93;^;^^ ^-^^-nm ^-^^-nM 13.64+^;^^ 4.22 2.02 1.38 0.82 



Errors are quoted at the 90% confidence level for 1 parameter of interest (Ax^ = 2.71); Fe narrow line energy fixed to E^.^ = 6.4 ke'V, width fixed to eV. ^ First line, intrinsic power 
law; second line, scattered power law. ^ Parameters for the states tied in the fit. " In units of 10"* photons keV~' cm"^ s^' @ 1 keV. ^ In units of 10^^ cm"^. ' Ionization parameter 
f = 4nF^/n (erg cm s"'), where Fjjj is the ionizing flux and n is the hydrogen nucleus density (partem"') of the illuminated slab. Complement to one of the covered-to-uncovered 
flux ratio. ' Line centroid energy, in keV. f Line width, in eV. ^ Ratio between the scattered flux and the partially absorbed component. 



components and different absorbers. Both models are statisti- 
cally acceptable. We then explored the possibility to explain the 
variability observed between November 2002 and January 2006 
(both in shape and intensity) by changing only the best fit pa- 
rameters (but assuming the same components in the model): we 
applied the models described in §5.2 to the November 2002 and 
(B) datasets. The two models provide a good description for the 
two spectra with reasonable values for the model parameters. In 
the absorption-based model, the soft gaussian emission fine at 
~ 0.57 keV is not requested by the November 2002 dataset: its 
EW upper limit is 273 eV at 90% confidence level. 

To conclude the analysis, we tested the possibiUty of fitting 
these models to the three states simultaneously, changing as less 
parameters as possible. 

The large number of variables involved in the reflection- 
based model (some of them degenerate in the fit) can make 
difficult the identification of the parameters responsible of the 
observed variations. To choose the parameters to be free, we 
started by fitting the spectra simultaneously with all but the nor- 
malization tied toghether. We then thawed one parameter at a 
time, searching for an improvement in the fit. The same step 
was repeated with an increasing number of free parameters, in 
various combinations, until the improvement stops to be signif- 
icant (i.e., F-test probability lower than 95%). In the final fit, 
the parameters identified in this way were kept free, while we 
tied the others between the states. Variations in the irmer ionized 
absorber, partially covering the central source, coupled with an 
increase of the reflection fraction, produce the differences ob- 
served between (A) and (B) spectra. To explain the change be- 
tween November 2002 and January 2006, both absorbing layers 
must vary; a lower accretion disk ionization state in the oldest 
observation is also required. 

The absorption-based model can reproduce the three spec- 
tra by varying three parameters: the intrinsic power law nor- 
malization, the covering fraction of the iimer absorber and the 
ionization parameter of the outer absorber. The power law nor- 
malization varies roughly proportionally to the observed hard 
(3 - 10 keV) flux. Variations observed between (A) and (B) 



datasets can be explained by means of a different covering 
fraction of the inner ionized absorber. The differences between 
November 2002 and January 2006 spectra can be attributed to a 
different covering fraction of the irmer absorber plus a variation 
of the ionization state of the outer absorber. 

The adopted models are able to describe the three states in 
a consistent way (x^/d.o.f. = 379.6/364 and 417.5/376; see 
Fig. 5): the best fit value for the parameters are reported in 
Table 5 and Table 6. 



7. Summary and Discussion 

We have presented a detailed spectral analysis of aU the EPIC 
XMM-Newton data available of PG 1535+547. The source has 
been first observed in November 2002, and then re-observed in 
January 2006; the latter data are presented here for the first time. 

7.1. X-ray emission 

Several features are clearly evident in the spectra of 
PG 1535+547 (see Fig. 2): a strong absorption between ~ 0.7 
and 2 keV, a spectral flattening at lower energies, a complex 
spectral shape in the iron band, and a drop in the emission at 
higher energies. Moreover, strong variations in shape and inten- 
sity are evident both on short and long timescales. The flux at 
energies < 1 keV increases of a factor ~ 10 in three years, then 
decreases of a factor ~ 2 in six days. Above 3 keV the variabiUty 
is significantly less important, with a change in flux of a factor 
~ 2 and ~ 1.3 in the two time intervals, respectively. 

The three broad band (0.4- 10 keV) EPIC spectra were anal- 
ysed in the framework of two different models: an "absorption- 
based" scenario and a "relativistic reflection-based" one. In both 
cases, a warm absorber is detected, producing the strong absorp- 
tion features observed in the 0.7 - 2 keV energy range. The spec- 
tral curvature observed at high energy can be reproduced either 
by a partial covering absorber or by a blurred relativistic reflec- 
tion component. In both cases, also a narrow gaussian emission 
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line, consistent with Ka fluorescence from almost neutral iron, 
is detected. 

To account for the observed spectra assuming an accretion 
disk reflection contribution, the continuum emission must be in- 
tercepted by two layers of matter. The innermost, ionized ab- 
sorber covers a high fraction (between ~ 0.7 and ~ 0.9) of the 
emission. The second layer, fully covering the source, is almost 
neutral ~ 3.3 erg cm s"'). As noted, this model suggest a re- 
flection produced in regions penetrating deeply towards the cen- 
tral blak hole < 2.50 rg). The disk emissivity power law 
index is close to 4, in agreement with theoretical expectations in 
case of a small height of the primary source {e.g., Martocchia 
et al. 2002). The reflection fraction, defined as the ratio between 
the flux in the ionized reflection component and the total flux. 



WA ■ IpcWA ■ (pi + d^k) + line] NH ■ Ipl + WA ■ [line + [kWA ■ (pl + Una)]} 




Fig. 5 Unfolded EPIC pn spectra of PG 1535+547 with the two adopted 
models fitted simultaneously to the three observations (see §6): the 
January 2006 (A) spectra (filled circles and red lines, upper panels) 
and (B) spectra (open triangles and gray lines, middle panels), and the 
November 2002 spectra (stars and sky-blue lines, lower panels). The 
solid lines correspond to the total best fit models. Over-plotted are the 
different components of each model: left panels, reflection from the disk 
(dashed line), power law (dash-dotted line), and narrow emission line 
(dotted line), with the thickest lines identifying the fraction covered by 
the partial covering warm absorbed ([1] and [2] respectively), while the 
thinnest lines identify the fraction uncovered, ([3] and [4] respectively); 
right panels, scattered component (dashed line, [5]), absorbed power 
law (dash-dotted line) and emission lines (dotted line), with the thick- 
est lines identifying the fraction covered by the partial covering warm 
absorber, [6], and the thinnest the fraction uncovered, [7]. 



is of the order of 50 - 60%, relatively large with respect to the 
20 - 30% expected in the case of an isotropic illumination. The 
large amount of reflection from the disk observed in all states 
requires some mechanism to enhance the reflection component 
with respect to the intrinsic continuum, resulting in a strongly 
suppressed intrinsic emission. Recent theoretical works, based 
on the "fight bending" eff'ect (Miniutti et al. 2003; Miniutti & 
Fabian 2004), or an infiomogeneous structure of tfie accretion 
disk due to violent clumping instabilities (Merloni et al. 2006), 
predict tfiis kind of effect. A similar spectral befiaviour fias 
been recently observed in tfie XMM-Newton spectra of Mrk 841 
(Petrucci et al. 2007). 

Tfie main origin of tfie observed variability results to be a 
cfiange in tfie absorbers properties (column densities and cov- 
ering fractions), coupled witfi a cfiange in tfie reflected compo- 
nent and tfie disk ionization state (see Table 5). Tfie coverage 
of tfie inner absorber decreases passing from November 2002 
to (B) and (A) spectra. A fiigfier disk reflection fraction and 
a fiigfier coverage of tfie inner absorber explain tfie spectral 
cfianges observed passing from (A) to (B) states. Tfie spectral 
index of tfie intrinsic continuum (F ~ 2) is consistent witfiin 
tfie errors witfi tfie value found for radio-quiet AGN (F ~ 1.9; 
Piconcelli et al. 2005; Nandra et al. 2007) and tfie mean value 
found for otfier NLS 1 s in tfie 2 - 1 ke V energy range (F ~ 2.19; 
Leigfily 1999). Tfie differences between November 2002 and 
January 2006 spectra imply stronger variations. Tfie primary 
continuum cfianges botfi in sfiape (flatter in tfie oldest observa- 
tion, F ~ 1.8) and in intensity, and tfie accretion disk ioniza- 
tion state increases. Tfie column densities of tfie two absorbers 
cfiange in an opposite way: in November 2002 tfie column den- 
sity of tfie inner absorber is fiigfier (A^h ~ 2 x lO^-' cm"^), wfiile 
tfie outer absorber sfiows a column density lower tfian in tfie 
January 2006 spectra (A^h ~ 5 x 10^' cm"^). 

Finally, a naiTow Fe Ka emssion line is included in tfie 
model, probably produced by tfie obscuring torus postulated 
in AGN Unified Models (but otfier possible sites are tfie outer 
part of tfie accretion disk, or tfie optical BLR). Recently, a nar- 
row emission line associated witfi cold iron was found to be al- 
most ubiquitous in a sample of local radio-quiet Seyfert galax- 
ies observed witfi XMM-Newton (Nandra et al. 2007). However, 
~ 70% of tfiese objects present furtfier complexities in tfie iron 
band. Tfieir spectral analysis sfiows tfiat a model based on broad 
emission from an accretion disk is tfie best way to take into ac- 
count tfiis complexity. Jimenez-Bailon et al. (2005) noted tfiat 
a naiTow Fe Ka emission line is required in ~ 53% of tfie 
PG AGN sample studied by Piconcelli et al. (2005). Tfie detec- 
tion of broad relativistic lines is statistically significant in less 
tfian 10% of tfieir sources (namely PG 0007+106, PG 0050+124 
and PG 1116+215). Interestingly, tfieir fiard X-ray luminos- 
ity falls in tfie lower end of tfie range covered by tfie sample 
(L2-wkeV < ergs s"'): in tfie disk reflection framework, tfie 
unabsorbed 2-10 keV luminosity observed from PG 1535+547 
cfianges from 9.7 x lO'*^ ergs s"' to 1.2 x 10"*^ ergs s"' and 
7.2 X lO'*^ ergs s"' for tfie (A), (B) and November 2002 spectra, 
respectively. From tfieir XMM-Newton data analysis of a large 
sample of AGN, Guainazzi et al. (2006) found tfiat relativisti- 
cally broadened Fe Ka lines are significantly more common in 
low luminosity (L2-iokeV ^ 10^"' ergs s"') AGN. 

In tfie absorption-based scenario several components are 
needed to properly reproduce tfie data: tfie primary power-law 
emission is partially covered by tfie warm absorber above men- 
tioned, and totally covered by a second ionized absorber. A soft 
scattered component is also present. In addition to tfie Fe Ka 
line, an O vii emission line (not required in tfie reflection-based 
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model) is also detected at £ ~ 0.57 keV. Finally, the whole emis- 
sion is absorbed by a rather low (~ 2 x 10^' cm"^) column den- 
sity of neutral gas. 

The observed variability can be accounted for by variations 
of the two warm absorber components. Between November 2002 
and January 2006 observations, the covering fraction of the inner 
absorber decreased significantly (passing from ~ 0.9 to ~ 0.55), 
leading to a different shape of the complex curvature seen at high 
energies. In the meanwhile the ionization state of the outer ab- 
sorber increased (from ~ 60 to ~ 90 erg cm s ') in a manner 
roughly proportional to the increase in unabsorbed luminosity. 
This may suggest that the density and location of the warm ab- 
sorber has not changed in this lapse of time. Instead, the differ- 
ences observed between (A) and (B) spectra can be explained by 
a variation of the covering fraction of the inner ionized absorber. 

The inner, partially covering warm absorber presents a lower 
(^ ~ 20 erg cm s"') ionization state than the fully covering warm 
absorber. In accretion disk winds theoretical models (e.g., Proga 
& Kallman 2004), two phase of ionized gas are predicted to co- 
exhist in the inner regions of accretion disks, where the wind 
is "launched". In particular, the wind itself is shielded from the 
central ionizing source by a higher ionization component. In this 
framework, it is possible that we are looking PG 1535+547 at 
the base of the wind, with the partial covering absorber corre- 
sponding to the "wind blobs" and the totally covering absorber 
corresponding to the highly ionized shielding phase. The latter 
warm absorber component may be responsible for the scattering 
of the primary emission into our line of sight. The very small 
scattering fraction (few parts per thousand) could thus be ex- 
plained by the low column density (~ 4 X lO'^'^ cm"'^) of the 
scatterer, coupled with a possible anisotropic distribution around 
the source. In this scenario, the emission line at ~ 6.4 keV can 
be produced by transmission through the high column density, 
partially covering absorber, while the line detected at ~ 0.6 keV 
could be associated to photoionized emission occurring in the 
low column density, totally covering absorber. The presence of 
unresolved photoionized emission lines would then naturally ex- 
plain the steepening of the scattered power law with respect to 
the intrinsic one. The primary power law emission is somewhat 
steep (F ~ 2.5), but consistent within the errors with the values 
foundbyLeighly (1999). Finally, the small (A^H ~ 2x10^' cm^^) 
neutral column density absorbing the whole emission and af- 
fecting only the data below 0.5 keV could be ascribed either to 
dust lanes in the host galaxy or to the presence of dust in the 
wind and/or scattering region. Overall, the complex spectra of 
PG 1535+547 seem thus fully compatible with the presence of a 
strong and complex ionized outflow along the line of sight. 

7.2. The broad-band SED of PG 1535+547 

Considering the strong variabihty observed in the X-ray band, it 
is interesting to look at the multiwavelength emission properties 
of PG 1535+547. 

In the UV band, Gallagher et al. (2001) noted a variation of 
30% between the lUE data of November 1982 and May 1984, 
and earlier lUE observation from June 1982 and September 
1982, that are found consistend with later HST spectra. The 
fluxes observed by the OM simultaneously to the new EPIC 
spectra presented here are in good agreement with the HST data. 

At long wavelengths, PG 1535+547 is the only radio-quiet 
AGN observed by Neugebauer & Matthews (1999) where a vari- 
ation in the near-infrared (NIR) emission from 1980 to 1998 is 
reported: in particular, the flux in 4 bands from 1 .27/im to 3.7/im 
rises of 0.5 mag between 1993 and 1996, going back to the pre- 



vious value between 1996 and 1998. The continuum emission 
observed in Seyfert galaxies at A > 2 is very likely domi- 
nated by thermal radiation from the inner edge of the dust torus 
(e.g., Kobayashi et al. 1993), heated by the radiation originated 
in the central active nucleus (Haas et al. 2003). How much the 
torus extends towards the black hole depends on the irradiating 
luminosity and on the dust composition, that determine the dis- 
tance at which the temperature is so high that the dust begins to 
sublimate. So, the dust evaporation radius can be estimated from 

the bolometric luminosity, /?evap - 0.06 V'^boi/(10''^ergs s"') pc 
(Netzer & Laor 1993; Hill et al. 1996). Assuming that the inner 
radius of the torus roughly corresponds to 7?evap> and adopting 
the bolometric luminosity estimated from the [O ni] line width 
by Zhang & Wang (2006), Lboi ~ 1-4 X 10"^^ ergs s we can es- 
timate the dimension of the region producing the observed NIR 
emission, /?torus ~ 0.07 pc. The minimum timescale for an ob- 
servable variability implied by /?toi-us is of the order of a year. 
Therefore, it is consistent with the possibility that the observed 
NIR variability is the reaction of the torus to the long timescale 
X-ray variable emission. 

In §2 we have mentioned the high optical continuum polar- 
ization showed by PG 1535+547, P = 2.5% (Berriman et al. 
1990). The optical polarization spectrum presented by Smith 
et al. (1997) exhibits a rise in the linear polarization P toward 
shorter wavelength (roughly proportional to /I"'), that contin- 
ues into the UV band down to 2000 A. A complex behaviour 
of the polarization position angle across the Ha feature is also 
observed. To explain the continuum and line polarizations, the 
authors suggest the presence of different scattering regions, with 
at least part of the scattering material intermixed with the line- 
emitting gas. Warm electrons responsible for the observed polar- 
ization in the emission-Une wings (i.e., the interaction with the 
warm absorber required by the X-ray observations) would be 
a possible alternative. Considering the two-component scatter- 
ing geometries recently proposed to account for the polarization 
properties of aU Seyfert galaxies simultaneously (Smith et al. 
2002, 2004), we can tiiink that we are observing PG 1535+547 
at a "privileged" angle. Although the available informations do 
not allow Smith et al. (2004) to confidently assign a polarization 
class ("polar-" or "equatorially-" scattered Seyfert) to this ob- 
ject, it is possible that our Une of sight passes very close to the 
absorbing torus. This seems to be a coherent scenario that take 
into account observations in different bands; further multiwave- 
length observations are required to test this hypothesis. 

In Fig. 6 we present the broad-band SED of PG 1535+547. 
A stronger high-energy variability than at low frequencies is 
clearly evident. The light-green dot-dashed fines corresponds to 
A = 3000 A and E = 2 keV in the source rest frame, i.e. where 
the low and high energy fluxes used by Brandt et al. (2000) to 
evaluate the are calculated. The change observed in the X- 
ray flux in front of a roughly constant UV flux impUes a variabil- 
ity in the X-ray-to-optical index even higher than previously re- 
ported: Gallagher et al. (2001) estimated for ROSAT and ASCA 
observations q-qx < -2.17 and qtox = -2.03. The observed in- 
dices measured from the new XMM-Newton spectra are consid- 
erably higher: we can compute the optical fluxes rescafing the 
UV HST spectrum to the OM data (that, we stress, are simul- 
taneous to the X-ray fluxes), obtaining aox = -1.66 and —1.65 
for the (A) and (B) spectra. This result implies a cfiange in the 
classification of PG 1535+547 as X-ray weak source. 

Many authors have investigated the relation between rest- 
frame UV and soft X-ray AGN emission, and its dependence 
with redshift and/or optical luminosity. Most studies have con- 
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Fig. 6 Spectral energy distribution for PG 1535+547. At low frequen- 
cies, open triangles are infrared and optical continuum data from 
Neugebauer et al. (1987), and the thin solid curve is the UV spectrum 
from HST (Gallagher et al. 2001). Yellow filled stars indicate the XMM- 
Newton OM fluxes observed in January 2006. In the X-ray energy range, 
ROSAT data point and upper limit are marked with open squares, while 
the black dotted line indicates the ASCA best fitting model (Gallagher 
et al. 2001). The three observed XMM-Newton EPIC fluxes are shown 
with the black long dashed line, January 2006 spectrum (B); red solid 
line, January 2006 spectrum (A); blue dashed line, November 2002 
spectrum. The vertical thin dot-dashed light-green lines correspond to 
/Iff = 3000 A and = 2 keV, while the green dot-dot-dot-dashed line 
corresponds to = 2500 A. 



eluded that there is no evidence for a redshift dependence, while 
the X-ray emission resuhs to be correlated with the UV emis- 
sion, and the ratio of the monochromatic X-ray and UV lumi- 
nosities (Tox decreases with increasing the latter (e.g., Vignali 
et al. 2003; Strateva et al. 2005; Steffen et al. 2006; Just et al. 
2007). In particular, Strateva et al. (2005) combined radio quiet 
sources from the SDSS, a heterogeneous low-redshift Seyfert 1 
sample, and a heterogeneous high-redshift sample. Performing 
their analysis both with and without the Seyfert 1 sample, 
they found that the monochromatic luminosities at 2500 A 
and 2 keV are correlated (with a slope lower than 1). The 
broad-band spectral index, calculated here at /} = 2500 A, 

"ox = log(f'2kev/^'2500A)/log(^2kev/v2500A)' results to be anti- 
correlated to the rest-frame monochromatic UV luminosity. The 
January 2006 XMM-Newton observations put PG 1535+547 
outside both relations (see Fig 7). From the observed luminos- 
ity at 2 keV we would expect a spectral index'"* a'^^ ~ -1.09, 
while the UV luminosity obtained from the OM data would im- 
ply an index a'^^ 1.29. Both values differ considerably from 

the observed indices, a'^^ - -1.66 and -1.67 for the (A) and 
(B) spectra, respectively. These numbers should be taken with 
caution: the host-galaxy contribution was not excluded from the 
January 2006 OM data, thus the nuclear UV flux could be lower 
than the value estimated here. However, taking this considera- 
tion in mind, PG 1535+547 appears to be still underluminous in 
the X-ray band if compared with optically selected sources in the 
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Fig. 7 Dependence on the 2500 A monochromatic luminosity of a'^^ 
(upper panel) and L2keV (lower panel). Black filled circles mark data 
from Strateva et al. (2005) for their "main" SDSS sample (155 objects 
with 0.1 < z < 4.5); arrows indicate upper limits in the X-ray detec- 
tion. Dashed lines are the best-fit linear relations for their combined 
sample ("main" SDSS sample + high-redshift sample + Seyfert 1 sam- 
ple): a'g^ = -0.136 ■ logLjgooA 2.616 (upper panel) and logL2keV = 
0.648 • log ^2500 A 6.734 (lower panel). Red open star marks the posi- 
tion of PG 1535+547 as observed in January 2006 (to avoid confusion 
we report only the values from the (A) data, being the two states indis- 
tinguishable). 

same bin of UV luminosity'^, a property shown by other NLSls 
with high-energy spectral complexity (e.g., Gallo 2006). 



8. Conclusions 

PG 1535+547 is a bright NaiTow Line Seyfert 1 galaxy, whose 
high energy emission shows strong variability both in shape 
and flux. In this paper, we have reported a consistent analy- 
sis of the three available XMM-Newton datasets of this object; 
the observations taken in January 2006 are presented here for 
the first time. The availability of high quality data obtained at 
different epochs with the same satellite opened us the very in- 
triguing possibility to study in details the complex emission of 
PG 1535+547, the origin of the spectral features as well as of 
the strong variability observed. The main aspects of this work 
can be summarized as follows. 

1 . Strong flux and spectral variability is observed on timescales 
of years and shorter; the strongest variation is observed in 
the soft X-ray band, at energies lower than 3 keV. The im- 
portance of absorption structures is evident in all the states, a 
result consistent with the strong polarization continuum and 
the absorption features observed in the optical-UV. 



■"2500 A 



1.65 X 10^' and 2 x 10 ergs s 'Hz respectively. 



'5 From Fig. 10 in Strateva et al. (2005), only 1 of 37 Seyfert 1 
falls in this bin of a'^^ (upper panel), while in their main SDSS sam- 
ple they have only 1 AGN with ^2500 A ^ ^-^^ ^ ^^S,^ ^ ' ' 
with -1.7 < Q-Qx < -1.6 (compared with 16 sources with i2500A ^ 
3.16 X 10'" ergs s"' Hz"', lower panel). 
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2. The broad band (0.4 - 10 keV) EPIC spectra can be well 
described either assuming an intrinsic continuum and reflec- 
tion from an accretion disk covered by a system of ionized 
absorbing materials (the "reflection-based" model), or by an 
intrinsic continuum and a scattered component, with a com- 
plex combination of neutral absorbing material and a two- 
phase warm absorber covering the central source in different 
ways (the "absorption-based" model). In the former case, 
the XMM-Newton observations imply evidence of impor- 
tant reprocessing from the disk of the intrinsic continuum 
(in terms of relativistically blurred reflection component and 
broad iron line) for aU the three states. 

3. In both scenarios, the analysis suggests that the observed 
variability can be ascribed mainly to changes in the ab- 
sorbers, although in the former the ionization of the emit- 
ting region of the disk plays a role, too. Further deeper 
observations are needed to model correctly the continuum 
in the 4-7 keV range and better constrain the parame- 
ters of the broad line and/or the absorbing structures, a task 
made particularly difficult by the well known interplay be- 
tween absorbers and broad features in the resulting spec- 
trum. Moreover, a spectral coverage extended at higher en- 
ergy would be able to fix with higher accuracy the intrinsic 
spectral index and the contribution from the reflection. 

4. Finally, unrelated X-ray and optical variabilities imply a 
change in the classification ofPG 1535+547 as X-ray weak 
source. The classification of a source as X-ray weak AGN 
seems to depend strongly on the date of the high energy ob- 
servations; the suggestion is that a similar behaviour can oc- 
cur also for other sources in intermediate range -2 < otox < 
-1.7. The nature of these sources as homogeneous class is 
not so clear, numbering among them objects like PHL 1811, 
that seems to be intrinsically X-ray weak (Leighly et al. 
2007a,b), or PG 2112-1-059, for which the weakness and the 
X-ray variability are ascribed to change in the intrinsic con- 
tinuum instead of to absorption effects (Schartel et al. 2007). 

To conclude, whatever the model used, the changes observed 
in the emission of PG 1535-1-547 imply physical variations oc- 
curring in regions close to the central black hole. Moreover, aox 
is definitely a variable parameter (see also Gallo 2006; Grupe 
et al. 2006): it seems that some sources go through phases where 
they appear to be X-ray weak and for most of them the X-ray 
weakness is a temporary event. Further investigations of their 
high energy emission with a simultaneous multiwavelength cov- 
erage are needed. 
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Fig. 4 Comparison of the unfolded best fits to the EPIC-pn January 2006 (A) data of PG 1535+547 with the two adopted models (see §5.2; upper 
panels). The residuals in terms of sigmas are plotted in the lower panels. The solid lines correspond to the total best fit models; over-plotted 
are the contributions of the different components. The unabsorbed emission, labeled as [Unabs.], is also shown. Left panels: Reflection from the 
disk (dashed line) and intrinsic power law (dash-dot-dotted line) partially covered by an warm absorber (the thickest lines identify the fraction 
of continuum covered, [1] and [2] respectively; the thinnest lines identify the fraction uncovered, [3] and [4] respectively), with a narrow line 
superimposed (dash-dotted line); the whole continuum was seen through a secon ionized layer of material. Right panels: Scattered component 
(dashed line, [5]), absorbed power law (dash-dot-dotted line) and narrow emission lines (dash-dotted line), partially covered by an warm absorber 
(the thickest line identify the fraction of continuum covered, [6], the thinnest the fraction uncovered, [7]); the whole continuum was seen through 
layer of neutral material. 



